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PHASE/GAIN IMBALANCE ESTIMATION OR COMPENSATION 
BACKGROUND 
Technical Field 

The present invention relates to a receiver for estimation or compensation of 
5 phase imbalance or gain imbalance and methods which can be applied in a receiver for 
estimation or compensation of phase imbalance or gain imbalance. A transmission 
technique used in the receiver is based on a QPSK (quadrature phase shift keying) 
modulation and a modulation scheme based on complex scrambling code. 

Description of the Related Art 
10 The Universal Mobile Telecommunications System or UMTS is one of the 

major new 'third generation' (3G) mobile communications systems being developed within 
the framework defined by the ITU (International Telecommunication Union) and known as 
IMT-2000. 

Methods for estimation or compensation of phase imbalance or gain 
15 imbalance typically rely on reference signal generation, adaptive interference cancellation, 
blind source separation or a least square approach in the Fourier domain. Such methods are 
known, e.g., from "Digital I/Q imbalance compensation in a low-if receiver", Jack P. Glas 
(Bell Labs)-1998-IEEE, "Improving QPSK Demodulator Performance for Quadrature 
Receiver with Information from Amplitude and Phase Imbalance Correction", Hung Ngyen 
20 (The Aerospace Corp.)-2000-IEEE, "Blind Source Separation Based I/Q imbalance 
compensation", Valkama, Renfors, Koivunen-2000-IEEE, "Advanced methods for I/Q 
imbalance compensation in Communication Receivers", M. Valkama, M. Renfros, V. 
Koivunen-IEEE transactions on signal processing- Vol. 49, no. 10, Oct. 2001. 

WO 0150616 discloses a receiver for correcting the phase error in a received 
25 signal. However, this receiver does not provide an efficient estimation or compensation of 
phase imbalance or gain imbalance in the received signal. 
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BRIEF SUMMARY 

It is an object of the present invention to provide an efficient estimation or 
compensation of the phase imbalance or the gain imbalance. 

According to an exemplary embodiment of the present invention, this object 
5 is solved with a receiver for estimation or compensation of phase imbalance or gain 
imbalance utilizing a QPSK modulation scheme and a modulation scheme based on 
complex scrambling code, the demodulator comprising means for estimating the phase 
imbalance or gain imbalance before synchronization. 

Such a receiver according to an exemplary embodiment of the present 
10 invention as set forth in claim 1, allows the estimation or compensation of phase imbalance 
or gain imbalance that can affect the overall performance of the receiver, at an early stage, 
i.e., before synchronization, by e.g., exploiting the properties of the complex scrambling 
code. Thus, phase imbalances or amplitude imbalances will not be able to introduce losses 
in the further phases of the connection. Furthermore, if the bandwidth inside which time 
15 variations in the imbalances occur is lower than the transmission bandwidth, the 
demodulator is also able to track these time variations. 

The phase/gain imbalance estimation and compensation according to the 
present invention is e.g., tailored to the UMTS modulation scheme, i.e., may exploit the 
special properties of the UMTS modulation scheme for estimation and compensation of 
20 imperfections like phase-gain imbalances in the RF receiver or demodulator in an early 
stage of the receiver. 

According to another exemplary embodiment of the present invention as set 
forth in claim 2, 1 and Q components of the received signal affected by phase imbalance or 
gain imbalance after demodulation are determined. A ratio between a cross correlation of 
25 the I and Q components (<IQ>) and the mean value of the square of the I component (<I 2 >) 
or a ratio between a cross correlation of the I and Q components and a square root of the 
product between the mean value of the square of the I component and the mean value of 
the square of the Q component ((<I 2 ><Q 2 >) 1/2 ) or a ratio between the mean value of the 
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square of the Q component, (<Q 2 >) and the mean value of the square of the I (<I 2 >) 
component are determined. 

According to yet another exemplary embodiment of the present invention as 
set forth in claim 3, a low pass for low pass filtering the signals are being implemented. By 
5 doing so, errors in the estimation of phase imbalance or gain imbalance can very efficiently 
be reduced, while computing mean values. 

According to yet another exemplary embodiment of the present invention as 
set forth in claim 4, the receiver comprises means for compensating the phase imbalance or 
gain imbalance before synchronization. 
10 According to another exemplary embodiment as set forth in claim 5, the 

receiver is a WCDMA (UMTS) receiver (Wide Band Code Division Multiple Access 
receiver). Thus, as a basic principle exploited in the WCDMA (UMTS) receiver, the 
channel is frequency selective with respect to the transmission bandwidth allowing for a 
very accurate estimation of the phase/gain imbalance before synchronization. 
15 According to another exemplary embodiment as set forth in claim 6, the 

estimation of the phase imbalance or gain imbalance is carried out iteratively providing for 
very accurate results. 

According to yet another exemplary embodiment as set forth in claim 7, a 
method for estimation or compensation of phase imbalance or gain imbalance in a receiver 
20 is provided, which may advantageously exploit the properties of the complex scrambling 
code and of the overall UMTS modulation scheme in the early stage of the receiver by 
estimating the phase imbalance or gain imbalance before synchronization. 

According to an exemplary embodiment of the present invention as set forth 
in claim 8, a very simple method is provided, requiring only minimized computation 
25 power, and which may advantageously allow an iterative solution for estimating 
compensating both phase and gain. 

According to yet another exemplary embodiment of the present invention as 
set forth in claim 9, a feed-forward scheme or a feed-back scheme is established. 
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According to another exemplary embodiment as set forth in claim 10, the 
estimation of the phase imbalance or gain imbalance is carried out iteratively providing for 
a very efficient compensation/estimation. 

According to the present invention, there is also provided a computer 
5 program for estimation or compensation of phase imbalance or gain imbalance as set forth 
in claim 1 1 . 

It may be seen as a gist of the present invention that the phase imbalance or 
gain imbalance is estimated or compensated before symbol synchronization. Thus the 
phase imbalance and gain imbalance will not be able to introduce losses in the further 
10 phases of the connection; also the method is able to track temporal drifts in phase and in 
gain if the bandwidth inside which these phenomena occur is lower than the transmission 
bandwidth. The advantage of this method is that it can be used both for phase imbalance 
estimation and compensation and gain imbalance estimation and compensation in the case 
that there is also a gain mismatch. 

1 5 BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

These and other aspects of the present invention will be apparent from and 
elucidated with reference to the embodiments described hereinafter. These embodiments 
will be described with reference to the following Figures: 

Figure 1 shows a feed-forward scheme for phase imbalance estimation and 
20 compensation according to an aspect of the present invention. 

Figure 2 shows a feed-back scheme for phase imbalance estimation and/or 
compensation according to an aspect of the present invention. 

Figure 3 shows a feed-back scheme for gain compensation according to an 
aspect of the present invention. 
25 Figure 4 shows schematically preferred application points for the method 

according to the present invention. 

Figure 5 shows simulation results for different filter lengths and a velocity 

of 3 km/h. 
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Figure 6 shows simulation results for different filter lengths and a velocity 

of 120km/h. 

Figure 7 shows simulation results for different filter lengths and a velocity 

of250km/h. 

5 Figure 8 shows the effect of the phase imbalance on the BER. 

Figure 9 shows simulation results for the square root error for 5 degrees 
initial phase error vs. low pass integration length. 

Figure 10 shows an estimation for a time varying phase (5 Hz). 
Figure 1 1 shows an estimated phase in accordance with an iterative 
10 implementation of the method according to the present invention. 

Figure 12 shows a gain error in accordance with an iterative implementation 
of the method according to the present invention. 



DETAILED DESCRIPTION 

In the following description of the above Figures, the same reference 
15 numbers are used for the same or corresponding elements. 

Figure 1 shows a device with a feed-forward scheme for phase imbalance 
estimation and compensation according to an aspect of the present invention. As will be 
shown below, the demodulation method according to which the device of Figure 1 is 
operating relies on the fact that the correlation matrix of the transmitted components, (i,q) 
20 is diagonal. As a consequence, if I and Q are the in-phase and quadrature components of 
the received signal (with I = k • i and Q = k • (-i-sincp + q-coscp) as shown in equations (1)), 
the correlation between I and Q, even after propagation through the fading multipath 
scenario, comprises only a term due to phase imbalance. This is due to the property of the 
QPSK modulation scheme and the complex scrambling code. This is shown in more detail 
25 in the following with respect to a UMTS system: 

The waveform transmitted in accordance with UMTS 
The UMTS transmission scheme in the downlink is based on a: 

- QPSK modulation 
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- Complex scrambling code 

The QPSK modulation scheme basically maps each two chips to be 
transmitted to a symbol of a quaternary constellation. QPSK symbols can be tough as 
complex numbers, and the application of the complex scrambling code to the QPSK 
5 symbols represented in terms of multiplication between complex numbers. In other words, 
to scramble the symbol (at the instant n) s(n)=(sl(n),s2(n)), the symbol for the scrambling 
code element c(n)=(cl(n),c2(n)) has to be multiplied in accordance with the following rule: 
s-c=(sl(n)-cl(n) - s2(n)-c2(n),sl(n)-c2(n)+s2(n)-cl(n)) (1) 
The two sequences composing the complex scrambling code, i.e., the sl(n) 
10 and s2(n) are selected in a family of Gold codes and have the following properties in first 
approximation: 

sl(n)*sl(n i + w) = 0, m±0 (2) 

sl(n)*s2(n + m) = 0, Vw 
where * refers to the correlation operation. 
15 It is important to note that above equations (2) hold true in first 

approximation, and that the correlation interval can be extended to the full length of the 
code, but also to sub-intervals. 

Indeed, if the components i and q of the transmitted signal given from the 
equation (1) are considered, the following equations (3) can be derived: 
20 i(n) = si (n)-cl (h) - s2{n)-c2{n) 

q{n) = sl(n)-c2(n) + s2(n)-c\(n) (3) 

i{n)*q{n) = <i{n)-q{n + m)>=0 

These relationships hold true for the maximal length of the Gold sequences, 
but also for intervals included if the included intervals are selected to have sufficient 
25 length. In the case of the UMTS, for the bandwidth involved, it is common understanding 
that the sufficient length can be extended to about 256 chips, i.e., the duration of a pilot 
symbol. 
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The in-phase and the quadrature components of the received signal are uncorrelated 

In the following, it is shown that after propagation through a multi-path 
fading environment, if there is no phase imbalance, the in-phase and the quadrature 
components of the received signal are still uncorrelated. 
5 At this point, i.e., before decoding, the signal is a mixture of voice, date and 

control information. The signal at the receiver can be written: 

K0=Z,ft- e ^- s ( r - T /)=Zi(0 

where the channel attenuation, delay and phase rotation is considered for each path. In 

what follows, I and Q are the in-phase and quadrature component of the received signal. I; 
10 and Qi are the received components on the i-path, and i,{t - t,) and q,(t - z,) represent the i- 

esima replica of the in-phase and quadrature components of the transmitted signals i(t) and 

q(t), due to the path /. r, is the delay of the transmission path. 

The I and Q components of n(t) can be expressed: 

Ii = gi[ii(t-Xi)-cos((pO + qi(t-Xi)sin(cpi)] 
15 Qi = gi[ii(t-Ti)-(-sin(cpi)) + qi(t-Ti)-cos(cpi)] 

from which it follows for the elements of the auto correlation matrix : 

<Ii 2 > = <gi 2 > • <ii 2 > • <cos(cp0 2 > + <gi 2 > • <qi 2 > • <sin(cpi) 2 >- 

2<g; 2 > • <i;qi> • <sin(cpi)cos(cpi)> 

<Qi 2 > = <gi 2 ><ii 2 > • <sin(cpi) 2 > + <g; 2 > • <q; 2 > • <cos((pi) 2 >- 
20 2<g; 2 > • <i;qi> • <sin((pi)cos((pi)> 

<Ii Qi> = <gi 2 > • <i; 2 > • <-sin(cpi) cos(cpi) + <g; 2 > ■ <q; 2 > ■ <sin(cpi) cos((pi)>- 

<gi 2 > • <i ; qi> • <sin((pi) 2 > - <g 2 > • <iiq ; > • <cos((pi) 2 > 

It has to be noted that <i ; 2 > = <q; 2 > = P/2 holds true for each I, with P being 

the transmitted power for each symbol. 
25 Also,<i;-qi> = 0 holds true for the orthogonal property of the in-phase and 

quadrature components of the transmitted signal. 
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Also, the terms gi, representing the channel attenuation for each path i, 
depend on the time. However, their dependence is slow with respect to the period onto 
which the mean operations are extended, and it can be assumed that they are uncorrelated 
with the temporal variations in the transmitted signal. As a consequence, the channel terms 
5 in the mean operations can be separated and written as follows: 

<gi 2 > = gi 2 (t) 

The same separation is performed for sin((pi) and cos(cpi) giving the phase 
information for each path. 

It follows that 
10 <Ii 2 >=<Qi 2 >= gi 2 (t)*P/2 

i.e., the received power is determined by the transmitted power and the channel gain. 

Furthermore, it follows that 

<IiQi>= 0 

i.e., the signal received on one path has in-phase and quadrature components that are 
15 orthogonal. 

Assuming a normalized power (per bit) during transmission (i.e., P/2=l), a 
consideration of the correlation between in-phase and quadrature components on the same 
path yields the following matrix: 

(/;-) (i,-<d) (g-( t ) o ) 

(q, i,) (a 2 ) J "I o g '(oJ 

20 where: 

s 2 (t)=Zs-(t) 

The correlation between the terms due to different paths 

In the following, the correlations between the terms due to different paths 

are considered. 
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The cross correlation between elements of different paths can be written as 

follows: 

<li lj> =< gi-gj> • <i;,ij> • <cos(cpi) cos((pj)> + <gig> • (qiaj> ■ <sin((pi) sin((pj>- 
<gi,gj> • <iicy><sin(cpj) cos((pi)> + <gigj> • <qiij><cos((pi) sin(cpj)> 
5 <Ii- Q> = <gigj> • <ii,ij> • <-sin(cpi) cos((pj)> + <gigj> ■ <q;qj> ■ <sin((pi) cos(cpj>- 
<gi-gj> • <ii-qj><sin((pi) sin(q)j)> + <grgj> ■ <qi ij><cos(cpi) ■ cos(q)j)> 

where, in the temporal mean, the terms due to the channel, g; and (pi, and the 
statistical operations involving these terms have been separated. Similar expressions can 
be derived for<QrQj> and <Qj Ii>. 
10 Since: 

<ii-ij> = <qrqj> = <qrij> = 0 

for i ^ j, all these terms are zero both in the short term and in the long term. 

Therefore it can be concluded that the correlation matrix between the 
components I and Q at the receiver is diagonal. Furthermore, elements on the principal 

15 diagonal depend on the time varying channel. However, it can be derived from equation 
(7) stated below, and from the fact that the channel is varying slowly with respect to the 
temporal extensions on which the temporal means for computing the correlations are taken, 
that this dependency does not influence the estimation of the gain/phase imbalance, since it 
is present at both the numerator and at the denominator of the final expression for sin(cp). 

20 As shown in Figure 1 , in a feed-forward scheme for phase imbalance 

estimation and compensation according to an aspect of the present invention, an incoming 
complex signal r(t) is received at a receiver station not shown in Figure 1 and is further 
processed such that the complex signal r(t) is provided to input connections 1 and 2 of the 
receiver or demodulator 3. After that, the complex signal r(t) is demodulated. Between the 

25 input connections 1 and 2, the signal Acos(a>t) represents the waveform used in the 

receiver to demodulate the incoming signal. Signal I on line 4 and signal Q on line 5 are 
called the I and Q components of the demodulated complex signal r(t) in the following. 
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The I signal or component is the "in-phase" component of the signal r(t) and the Q signal 
or component is the quadrature component. 

The I and Q components affected by phase imbalance or gain imbalance and 
are provided to lines 6 and 7, which are connected to lines 4 and 5. The I and Q 
5 components are fed into means 8 for taking the cross correlation of I and Q, such as a 

multiplier. Additionally, the I component is fed into means 9 for taking the quadrature of I, 
such as a multiplier which multiplies I with I. After that, the resulting signals are taken 
through low pass filters 12, 13 by lines 10, 11, respectively. For estimation of the phase 
imbalance, a ratio between the cross correlation of I and Q and the mean value of the 

10 square of I is generated by means for generating a ratio 14, such as a divider, and 

consecutively passed through low pass filter 15. After that, the resulting signal is provided 
to a phase/gain compensator 17 via a line 16. The phase/gain compensator 17 compensates 
the phase imbalance or the gain imbalance of the I and Q components. Thereupon, the 
estimated and compensated I and Q components i c and q c , are taken by lines 18 and 19, 

15 respectively, to be fed into means for synchronization 20, such as a UMTS synchronizer. 
After synchronization, the signal is ready for further processing, e.g., in reference to the 
standard architecture for the UMTS receiver, it can be fed into the rake receiver to exploit 
time diversity and then to the channel decoder to correct/reveal residual errors. 

Reference character 21 designates a multiplier or squaring device connected 

20 to line 5. An output of the multiplier 21 is connected to a low pass 22. An output of the 
low pass 22 is connected to another divider 23. The divider 23 receives the output signals 
from the low pass 13 and the low pass 22 and outputs a ratio of the two input signals. The 
output of the divider 23 is connected to another low pass 24, which in turn is connected to 
the phase/gain compensator 17. 

25 The operating principle of the receiver of Figure 1 and a method for 

operating the receiver of Figure 1 according to an exemplary embodiment of the present 
invention will be described in the following. 
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Phase imbalance estimation 

The method for phase imbalance estimation according to an exemplary 
embodiment is based on the observation that after the demodulator, the I and Q 
components affected by imbalance can be written as function of the i and q components 
5 without imbalance as shown in equation (4): 
I = k-i 

Q = k-(-i-sincp + q-coscp) (4) 

Equations (4) hold true if the phase of the I-branch of the demodulator is 
aligned with the phase of the i component of the received signal. 
10 (i,q) is assumed to be the transmitted signal for UMTS, i.e., the signal 

without including the effects due to the propagation through the multi-path fading channel. 
From the complex scrambling code properties and from the modulation scheme adopted the 
following equations (5) can be inferred. 
<i-q> = <q-i> = 0 

15 <H> = <q-q> = P/2 (5) 
Thus the correlation matrix of (i,q) is diagonal as already described above 

and the terms on the principal diagonal depend from the transmission power. 
From (4) and (5) it can be written: 

<IQ> = -k-<i 2 >sincp 

20 <I 2 > = k<i 2 > (6) 
<IQ>is the output signal of the low pass 12 and<I 2 >is the output signal of 
the low pass 13. Both signals are input to the means for generating a ratio 14. 

As shown above, it can be inferred that: 
sincp = -<I Q> / <I 2 > (7) 
25 where<IQ>is defined as cross correlation of I and Q, <I 2 >as quadrature of I, 

and<Q 2 > as quadrature or square of Q. 

When propagating through the multi-fading environment, the ideal QPSK 
(quadrature phase shift keying) constellation at the transmitter is rotated/attenuated/delayed 
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for each path. This has an influence on the cross correlation properties of (i,q) before the 
receiver. Furthermore, the signal becomes intrinsically not stationary and the statistics of 
the (i,q) samples get more complicated by the fact that there can be correlations between 
delayed versions of the same signal. It has been shown above that these effects have no 
5 influence and the residual presence in the cross correlation in the received components (i r , 
q r ) is only due to the phase imbalance. 

For compensation of the phase imbalance, the resulting signal is put into 
phase/gain compensator 17, which applies the following equations: 

ic =1 

10 q c =(Q + I-sincp) / coscp (8) 
where: 

cos(p = ((</ 2 > 2 -<I Q> 2 )/<I 2 > 2 y (9) 
i c and q c are the output signals of the phase/gain compensator 17 on the lines 

18 and 19. 



15 Gain imbalance estimation in the presence of phase imbalance 

In the following, the method is extended to the case in which there is also an 

amplitude mismatch besides the phase imbalance. 

In case there is also an amplitude mismatch, the above equations can be 

written as follows: 
20 I = g(t)(k r i) 

Q = g(t)(k 2 -(-i-sin(p + q-coscp) 

<I,Q> = ki • k 2 • <g(t) 2 ><q 2 >sin((p) 

<IQ>= -h ■ k 2 <g(tf><i 2 >sin((p)=-ki ■ k 2 <g(tf><q 2 >sm(f) 

<Q- Q>=k 2 2 • <g(t) 2 >(<i 2 >cos(cp) 2 +<q 2 >sin((p) 2 )=k 2 ■ <q 2 > 
25 <I- I>=ki 2 <g(t) 2 ><i 2 > 

From these equations it follows: 
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<I,Q> / (<I 2 xQ 2 >) 1/2 = -sin(cp) 

<Q 2 >/<I 2 > = k 2 2 /ki 2 (10) 
Where the ratio<Q 2 >/<I 2 > is the output signal of the divider 23. After a low 
pass filtering by means of low pass 23, the ratio <Q 2 >/<I 2 > is provided to the phase/gain 
5 compensator 17. 

Note that the temporal coefficient g(t) representing the amplitude of the 
channel dependence on time can be simplified, because it appears both in the numerator 
and in the denominator of equation (10), which can be considered a generalization of 
equation (7). The phase factors are removed also, because they would anyway disappear in 
10 the auto-correlation matrix. 

Feed-back scheme for phase imbalance 

Figure 2 shows a feed-back phase imbalance estimation and compensation 
device according to an aspect of the present invention. As shown in Figure 2, an incoming 
complex signal r(t) is received at a receiver station not shown in Figure 2 and is further 

15 processed such that the complex signal r(t) is provided to input connections 1 and 2 of the 
demodulator 3. After that, the complex signal r(t) is demodulated. Between the input 
connections 1 and 2, A-cos(oot) represents the waveform used in the receiver to demodulate 
the incoming signal. 

The I component on line 4 and the Q component on line 5 of the 

20 demodulated complex signal r(t), which are affected by phase imbalance or gain imbalance, 
are put through phase/gain compensator 17 and afterwards provided to lines 18 and 19. 
The output signals of the phase/gain compensator 17 are designated with i c and q c The i c 
and q c signals are fed into means 8 for taking the cross correlation of the I and Q 
components, such as a multiplier by means of lines 6 and 7. Additionally, the i c signal is 

25 fed into means 9 for taking the quadrature or square of the I component, such as a 

multiplier which multiplies the i c with itself. After that, the resulting signals are provided 
to low pass filters 12, 13 via lines 10, 11, respectively. 
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For estimation of the phase imbalance, a ratio between the cross correlation 
of I and Q and the quadrature of I is generated by means for generating a ratio 14, such as a 
divider. After that, the resulting signal is integrated by means of an integrator 25 provided 
to line 16 and fed back into the phase/gain compensator 17 located between lines 4, 5 and 
5 18, 19, for compensating/estimating the phase imbalance or the gain imbalance of the I and 
Q components. Thereupon, the estimated and compensated I and Q components i c and q c , 
respectively, are fed back to lines 18, 19 to be fed into means for synchronization 20 such 
as a UMTS synchronizer. 

In the following, an exemplary embodiment of a method of the present 
10 invention for operating the device of Figure 2 and its operating principle will be described. 

In the feed-back phase imbalance estimation and compensation device 
according to the present invention and the method for operating the feed-back phase 
imbalance estimation and compensation device according to the present invention, the 
compensated output of the phase/gain compensator 17 feeds the imbalance estimation. 
15 Some simplifications have been introduced in the formulas for the compensation and for 
the phase estimation. This is possible because in the next iteration in the loop the phase 
error gets smaller and the approximation of the sinus with the value of the angle and the 
cosines with 1 is appropriate. The principle of the method will be described and explained 
in the following: 

20 The following equation can be used for the estimation of the phase cp: 



The phase compensation carried out in the phase/gain compensator 17 can 
be described with the following equation: 



cp = -<I Q> / <I 2 > 



(11) 



(12) 



25 q c =(Q + I-q>) 



(13) 



A further simplification which can be introduced for the phase 
compensation feedback method is to use the following expression for the phase: 
cp = -k*<TQ> (un-normalized expression for the error) 
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where k is a constant. This is possible because this term is equal to zero when the phase 
imbalance is zero. Hence, this the signal represented by this term can be then used as error 
in the feedback scheme. 

According to an aspect of the present invention, a more efficient way to 
5 implement the feedback scheme is via an iterative procedure. According to this aspect, the 
low pass filters and the integral may unified in a single computation step for 
implementation. This will be described with more detail in the following: 

Iterative implementation for phase imbalance 

In the following, i c (n), qdn) are the phase compensated received in-phase 
10 components and quadrature components. 

1) At an instant n, the following error function is computed: 

e =i c (n)-q c (n) 

2) Then, the following equation is used to obtain a filtered integrated expression for 
the error: 

15 e/=e + p-e 

where the parameter p is chosen by matching convergence speed and stability. 

3) Then, the next (subsequent) i c (n +1), q c (n +1) is computed from the next 
i{n +1), q(n +1) received samples as follows: 

i c (n +1), i(n +1) 
20 q c {n+\) = q(n+\) + i(n+\)-e f 

Then, the method returns iteratively to step 1). It has to be noted that the 
final value for ej represent the estimated sinus of the phase imbalance multiplied for the 
gain factor k. 

Feedback scheme for gain imbalance 
25 Figure 3 shows a device for gain compensation with a feedback scheme. 
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The signal i c on the line 18 is provided to means 30 for taking the quadrature 
or square via line 7. The output signal of the means 30 for taking the quadrature is filtered 
by means of a low pass 3 1 and then provided to means for subtracting 32 such as a 
subtractor. This signal is designated as signal A in Figure 3. 
5 The signal q c on the line 19 is provided to means 9 for taking the quadrature 

or square via line 6. The output signal of the means 13 for taking the quadrature is filtered 
by means of a low pass 3 1 and then provided to the means for subtracting 32. This signal 
is designated as signal B in Figure 3. The means for subtracting 32 are configured such 
that the output signal of the low pass 3 1 is subtracted from the output signal of the low pass 
1 0 which is indicated by B - A in Figure 3 . 

The device depicted in Figure 3 operates as follows: The received input is 
demodulated, and then gain compensated in the each iteration of the loop. The gain 
compensated components are squared, filtered and subtracted to build an error function 
which is integrated and then provided to the gain compensation block 17, to be processed 
15 together with the incoming signal. 

The equations that give a rationale of the method are the folio wings. 
Starting from the equations: 

<i-Q> ■ { \ 

I :— ^ = = sm(<p) 

V< I 2 x Q 2 > 

<Q 2 > = K> 
<I 2 > k 2 

20 wherein x can be defined as follows: 

x _ k 2 -k 2 _ <I 2 >-<Q 2 > „ 2-k x -{k x -k 2 ) _ 2-{k l -k 2 ) 
k 2 <I 2 > k 2 k x 

which becomes: 

x _ k 2 -k 2 _ <I 2 >-<Q 2 > ^ 2-k x \k x -k 2 ) _ 2-{k l -k 2 ) 
k 2 <I 2 > k 2 k x 

and 
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^= l -°- 5 - x =t 2 



(16) 



where the following approximation has been used: 



x = k l +k 2 = 



2-k, 



The gain compensated components are obtained as follows: 



5 i c = I 



q c = Q-g 



(17) 



Note that the derivation is based on an approximation becoming more 
precise with a decreasing difference between kj and under the condition that it is realized 
in the next iterations of the loop. 



x 1 = k\ - k\ (un-normalized expression for the error) 

for applying the method. In fact, this last expression is zero for gain balanced components, 
and does not involve the ratio operation which can present some problems in the presence 
of noise. 

15 Iterative implementation for gain imbalance estimation and compensation 



integration operation is represented by the following steps in accordance with an exemplary 
embodiment of the present invention: 

1) In an instant n, the following term is computed: 
20 e = i c (nf-q c (nf 

2) Then, the following equation is used to obtain a filtered integrated expression for 
the error: 

ej = e + p • e 

where the parameter p is chosen by matching convergence speed and stability 
25 3) Then, the gain g cab be derived as follows (k is a given constant): 



10 



A further simplification is obtained by considering: 



An even more efficient way to perform the low pass filtering and the 
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g = 1 - k • e f 

4) The computed gain is then applied to the next i and q samples: 



i c (n+\) = (n + \) 

q c {n+\) = q{n + \)-g 
5 5) Then, the method iteratively restarts from step 1) 

Note the final value of g will represent the estimated value for the gain. 
Advantageously, this iterative method is valid in presence of a phase imbalance and in 
absence of phase imbalance. 

Phase and gain imbalance: feedback/iterative method. 

10 The two loops represented in Figures 2 and 3 can be serially concatenated to 

obtain simultaneous compensation of phase and gain. In particular, the output of the loop 
represented in Figure 2 can be the input of the loop represented in Figure 3 and the final 
output would be compensated both in phase and in gain. This is a consequence of the fact 
that in the equations (10) phase and gain imbalances are independent. Also, if this 

15 independence does not apply in the case in which un-normalized expressions are used for 
the errors to reduce complexity, it can still be shown that the serial concatenation of the 
loops operates very efficiently. In fact, in this last case, the residual phase imbalance 
influence the estimation of the gain imbalance in the next iteration of the loops, but this 
influence becomes zero when the phase error becomes zero. Accordingly, the coupling 

20 created by using the un-normalized forms for the errors creates a reciprocal influence in 
terms of convergence speed, but does not affect the convergence itself. 

Figures 1 1 and 12 depict results achieved with two concatenated iterative 
loops in which the phase and the gain are simultaneously compensated. 

As will be shown in the following, the above iterative approach based on 

25 un-normalized errors allows for a very efficient and more general estimation and/or 
compensation of gain/phase imbalances. In fact, starting from the following equations: 
I = k r i (18) 
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Q = -k 3 i- sirup + k 2 q ■ cos(p 

in which £3 is used for the part of i injected on the q branch, it can be shown that: 
<I Q>= k x k 3 sin((p) < z 2 > 

< I 2 > - < Q 2 >=< i 2 > (k 2 -k 2 sen 2 <p)-klq 2 cos<p (19) 

5 As can be seen from these equations ( 1 9), the first one holds true if one of 

the two phase and gain imbalanced components is multiplied with an arbitrary factor. 
Also,<IQ>gives an error that is zero only if the phase imbalance is zero. On the other 
hand, in absence of phase imbalance, the second equation is zero, if and only if: 

< I 2 > - < Q 2 >=(< i 2 > k 2 - < q 2 > k 2 )=0; (20) 

io 

< q > kl 

In this way, a cascade scheme can be designed, in which the equation for the 
phase is applied on phase unbalanced components as given in equation (11). The output of 
this first loop feeds a second loop, which applies the procedure given in equations (12) and 
(13) to estimate and compensate for the gain imbalance. This scheme will compensate 

15 both phase and gain imbalance in a computational effective way whenever they can be 
expressed in the form given by the equations (18). 

Figure 4 shows components of a UMTS receiver. There is a radio frequency 
(rf) demodulator 40 which output signal is input to an A/D converter 41 . The output of 
A/D converter 41 is input to an RRC filter 42 which is connected to a down sampler 43. 

20 From the down sampler 43, the signal is output to a further processing. The lines 44, 45 
and 46 indicate positions within the UMTS receiver, where the above methods for phase 
and/or gain imbalance estimation and/or compensation may be applied. As can be seen 
from Figure 4, phase and/or gain imbalance estimation and/or compensation may be 
applied between the A/D converter 41 and the RRC filter 42, i.e., before the RRC filter 42, 

25 between the RRC filter 42 and the down sampler 43, i.e. , before the down sampler 43, and 
after the down sampler 43 before the further processing. Therefore, according to the 



19 



present invention, a great degree of freedom is provided for the design of a receiver 
according to the present invention. 

Performance of the inventive method and device 

Figures 5-7 show a ratio between the estimation error and the true value for 
5 filters of different length as resulting from a simulation based on typical conditions for 
UMTS receivers extracted from the standard specifications. The phase imbalance error 
reaching from 0 to 5 degrees and the length of the low pass 2 in Figure 1 are represented in 
the X, Y axes. The unit for the length of the low pass filter is 1 pilot symbol, i.e., 256 
chips. The length of the low pass filter 1 in Figure 1 is always assumed to be 256 chips. 
10 The Z-value E represents the square root of the variance of the estimated value of the actual 
imbalance. It is related to the mean magnitude of the residual phase error after 
compensation. 

In Figure 5, the conditions considered correspond to 4 paths, a velocity of 3 
km/h and lorx/loc =-3 dB. 
15 In Figure 6, the conditions considered correspond to 4 paths, a velocity of 

120 km/h and lorx/loc =-3 dB. 

In Figure 7, the conditions considered correspond to 4 paths, a velocity of 
250 km/h and lorx/loc =-3 dB. 

All three simulations show very good results for E for a filter length of low 
20 pass filter 2 in Figure 1 of more than 50 pilot symbols or a phase error of more than 2 
degrees. 

Figure 8 shows the effect of the phase imbalance on the BER (Bit-Error 
Rate) in a case in which AWGN propagation conditions (Additive White Gaussian Noise) 
are considered and no other losses are included. Between 0 and 4 degrees the BER 
25 increases rather slowly. From 4 to 7 degrees, the increase of the BER is moderate and from 
7 to 10 degrees, the BER increases fast. As can be seen in the range 0-5 degrees, the phase 
imbalance can imply an increase in the BER of about 60%. 
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Figure 9 shows the performance which can be achieved in typical 
conditions. The phase error decreases exponentially with the integration length. By 
choosing an integration length of 16 x 256 chips, corresponding to 16 pilot symbols, the 
phase error is reduced in mean of the 90%. 
5 Figure 10 shows results related to a case in which the phase is changing at a 

fixed frequency of 5 Hz. This case is used to give a representation of the tracking 
capability of the method. The conditions considered correspond to 4 paths and a velocity 
of 120 km/h. The effect of the noise is included, inter-cells interference at + 3 dB with 
respect to the desired signal. The integration period is 16 slots and the mean square error 

10 achieved is 0.57 on a peak phase error of 5 degrees. 

Figures 1 1 and 12 show the performances of the iterative procedure applied 
when both phase and gain imbalance arc present. Figure 1 1 shows an estimated phase with 
speed =3 km/h and 3 equivalent paths in the iterative application of the method at the 
location indicated with arrow 44 in Figure 4. The initial condition for the gain imbalance is 

15 0.8 and 5 degrees for the phase. The term ki in the equations (1 8) is equal to 0.5. 

Figure 1 shows the error for the gain with speed =3 km/h and 3 equivalent 
paths in the iterative application of the method at the location indicated with arrow 44 in 
Figure 4. The initial condition for the gain imbalance is 0.8 and 5 degrees for the phase. 
The term £ 3 in the equations (18) is equal to 0.5 

20 As can be taken from Figures 1 1 and 12, after approximately 300 pilot 

symbols, i.e., 2 frames corresponding to an interval of 20 msec, errors of 5 degrees in phase 
and 0.8 in gain are reduced to about 0.02 degrees for the phase and 0.01 for the gain. 

The method and device of the present invention is particularly well suited 
for UMTS/WCDMA, where a complex scrambling code is used to modulate the 

25 information symbols. However, the method and device of the present invention may also 
be applied for any transmission with QPSK modulation, a modulation scheme based on a 
complex scrambling code and a frequency selectivity of the channel with respect to 
transmission bandwidth. 
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The device according to the present invention may be designed by means of 
EPLDs, such as the Altera EPLDs or by means of suitable calculation devices. The 
methods of the present invention may be implemented as software program code in a 
suitable programming language such as C++ and stored on a suitable computer readable 
5 storage device such as a compact disc. 

The various embodiments described above can be combined to provide 
further embodiments. All of the U.S. patents, U.S. patent application publications, U.S. 
patent applications, foreign patents, foreign patent applications and non-patent publications 
referred to in this specification and/or listed in the Application Data Sheet, including but 
10 not limited to [insert list], are incorporated herein by reference, in their entirety. Aspects of 
the embodiments can be modified, if necessary to employ concepts of the various patents, 
applications and publications to provide yet further embodiments. 

These and other changes can be made to the embodiments in light of the 
above-detailed description. In general, in the following claims, the terms used should not 
15 be construed to limit the claims to the specific embodiments disclosed in the specification 
and the claims, but should be construed to include all possible embodiments along with the 
full scope of equivalents to which such claims are entitled. Accordingly, the claims are not 
limited by the disclosure. 
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